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GHOSTS OF BUMP ATTRACTORS IN STOCHASTIC NEURAL 
FIELDS: BOTTLENECKS AND EXTINCTION 


ZACHARY P. KILPATRICK* 

Abstract. We study the effects of additive noise on stationary bump solutions to spatially 
extended neural fields near a saddle-node bifurcation. The integral terms of these evolution equations 
have a weight kernel describing synaptic interactions between neurons at different locations of the 
network. Excited regions of the neural field correspond to parts of the domain whose fraction of active 
neurons exceeds a sharp threshold of a firing rate nonlinearity. For sufficiently low firing threshold, a 
stable bump coexists with an unstable bump and a homogeneous quiescent state. As the threshold is 
increased, the stable and unstable branch of bump solutions annihilate in a saddle node bifurcation. 
Near this criticality, we derive a quadratic amplitude equation that describes the slow evolution of 
the even mode (bump contractions) as it depends on the distance from the bifurcation. Beyond the 
bifurcation, bumps eventually become extinct, and the time it takes for this to occur increases for 
systems nearer the bifurcation. When noise is incorporated, a stochastic amplitude equation for the 
even mode can be derived, which can be analyzed to reveal bump extinction time both below and 
above the saddle-node. 

Key words. Stochastic partial differential equations, Langevin equation, Perturbation theory. 
Amplitude qquations, Saddle-node bifurcation 


1. Introduction. Continuum neural fields are a well-accepted model of spa- 
tiotemporal neuronal activity evolving within in vitro and in vivo brain tissue [6,11]- 
Wilson and Cowan initially introduced these nonlocal integrodifferential equations 
to model activity of neuronal populations in terms of mean firing rates [54]. While 
they discount the intricate dynamics of neuronal spiking, these models can qualita¬ 
tively capture a wide range of phenomena such as propagating activity waves observed 
in disinhibited slice preparations [27,44,45,47]. Neural field models exhibit a wide 
variety of spatiotemporal dynamics including traveling waves, Turing patterns, sta¬ 
tionary pulses, breathers, and spiral waves [14,17,19,31,38]. A distinct advantage of 
utilizing these continuum equations to model large-scale neural activity is that many 
analytical methods for studying their behavior can be adapted from nonlinear partial 
differential equations (PDEs) [6] . Recently, several authors have explored the impact 
of stochasticity on spatiotemporal patterns in neural fields [8, 30, 35] by employing 
techniques originally used to study stochastic front propagation in reaction-diffusion 
systems [48]. Typically, the approach is to perturb about a linearly stable solution 
of the deterministic system, under the assumption of weak noise. However, some re¬ 
cent efforts have been aimed at understanding the impact of noise on patterns near 
bifurcations [30,36]. 

In this work, we are particularly interested in how noise interacts with stationary 
pulse (bump) solutions near a saddle-node bifurcation at which a branch of stable 
bumps and a branch of unstable bumps annihilate [1] . Bumps are commonly utilized 
as a model of persistent and tuned neural activity underlying spatial working memory 
[21,55]. This activity tends to last for a few seconds, after which it is extinguished, 
to allow for subsequent memories to be formed [23]. One possible way to terminate 
these sustained activity patterns is by transiently synchronizing the spiking patterns 
of excitatory neurons that participate in the signal [25]. Another proposed mechanism 
for terminating persistent activity is a strong and brief global inhibitory signal, which 
would drive the system from the stable bump state to a stable uniform quiescent 
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state [10]. In terms of neural field and spiking models, this can be thought of as 
momentarily raising the firing threshold of the system, temporarily driving it beyond 
the saddle-node bifurcation from which the stable bump emerges. 

We focus on a scalar neural field model that supports stationary bump solutions 
for appropriate choices of parameters and constituent functions [ 1 , 11 ]: 

= -u{x, *) + - y)fiu{y, i))dy (1.1) 

where u{x, t) is the total synaptic input arriving to location x and time t, and w{x — y) 
describes the strength (amplitude) and polarity (sign) of synaptic connections from 
neurons at location y to neurons at location x. We assume w{x) is an even-symmetric 
function w{x) = w{—x) with a bounded integral f^w(x)dx over the spatial domain 
X € = {—Xoo, Xoo)- The nonlinearity f(u) is a firing rate function, which we take to 

be the sigmoid [54] 

^ 1 _|_ Q-v{u-e) ’ 

and we also find it useful to take the high gain limit 77 —^ oo, in which case: 

/(")="(“-»)={ i (13) 

allowing for analytical tractability in several of our calculations. It is important to note 
that ( 1 . 1 ) neglects several known features of neuronal networks including spike rate 
adaptation [26], propagation delays [29], synaptic depression [33], and refractoriness 
[16]. Thus, we assume we are focusing on a network where these effects are weak 
enough as to not impact our main results. 

Amari was the first to analyze (1.1) in detail, showing that when f{u) is defined 
to be a Heaviside function (1.3), the network supports stable stationary bump solu¬ 
tions when the weight function w{x) is a lateral inhibitory (Mexican hat) distribution 
satisfying: (i) w{x) > 0 for cc S [ 0 ,a;o) with w(xo) = 0 ; (ii) w(x) < 0 for a; G (a;o,a:oo); 
(iii) w(x) is decreasing on [ 0 ,a:o]; and (iv) w(x) has a unique minimum on [ 0 ,a;oo) at 
X = Xi with xi > xo and w(x) strictly increasing on (xi,Xca) [!]• Based on restrictions 
(i)-(iv), Amari made use of the integral of the weight function 

W{x)= [ w{y)dy (1.4) 

Jo 

to prove some of the main results of his seminal work. For instance, it is clear that 
kF(0) = 0 and W{x) = —W{—x) based on the above assumptions. Moreover, there 
will be a single maximum of the function W{x) on the interval (0,a:oo) given at 
X = xq, i.e. Wmax = Ta&yixW {x) = IF(a;o), due to conditions (i) and (ii), and 
w(xo) = 0. When 0 < W{xo) there are two bump solutions: one stable and one 
unstable (up to translation symmetry), and when 6 > W{xq) there are no bump 
solutions to (1.1). When 9 = 9(. = W{xo), there is a single marginally stable bump 
solution. It is at this point that the two branches (stable and unstable) of bump 
solutions meet and annihilate in a saddle-node bifurcation (Fig. 2.1). Dynamics of 
(1.1) for values of 9 beyond this saddle-node bifurcation evolve to quasi-stationary 
solutions resembling the ghost of the bump at 9^, lasting for a period of time inversely 
related to \/\9J^l)^\ [51]. A principled exploration of these dynamics (section 2) is 
one of the primary goals of this paper. 
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Fig. 2.1. Saddle-node bifurcation of bumps in (1.1) with a Heaviside firing rate function (1.3). 
(A.) Difference of Gaussians weight function w{x) = e~^ — has a Mexican hat profile with 

A = 0.4 < 1 and cr = 2 > 1. The critical bump half-width ac at the saddle-node satisfies the relation 
w(2ac) = 0. (^) The weight function integral (1.4) determines the bump half-widths a. When 9 is 
below the critical threshold 9c at the saddle-node, there are two stationary bump solutions to (1.1): 
one stable as and one unstable Ou- When 9 > 9c, there are zero equilibria, but the dynamics of (1.1) 
are slow in the bottleneck near Uc{x). 

As mentioned, the neural field equation (1.1) in the absence of noise has been 
analyzed extensively [1,11,17]. We expand upon these previous studies by also ex¬ 
ploring the impact of noise on stationary bump solutions to (1.1) near a saddle-node 
bifurcations (section 3). Additive noise is incorporated, so that the evolution of the 
neural field is now described by the spatially extended Langevin equation [4,6,30,40]: 


(iu{x, t) 


—u(x, t) 



y)f{u{y,t))dy 


dt -I- eAW{x, t), 


(1.5) 


where the term dW(a;,t) is the increment of a spatially varying Wiener process with 
mean defined by {dW{x,t)) = 0 and correlations dW{x,t)dW{y, s)) = C{x — y)d{t — 
s)dtds and e describes the amplitude of the noise, assumed to be weak (e <C 1). The 
function C{x — y) describes the spatial correlation in each noise increment between 
two points x,y G fl. 

2. Slow bump extinction in the deterministic system. We begin by ex¬ 
amining the dynamics of stationary bump solutions near a saddle-node bifurcation, 
where a stable and unstable branch of solutions annihilate. Our initial analysis focuses 
on the noise-free case W (x, t) = 0, allowing us to derive an amplitude equation that 
approximates the evolution of the bump height. Linearization of bumps in (1.1) typi¬ 
cally reveals that they are marginally stable to translating perturbations, so the overall 
stability is characterized by the stability to even perturbations that expand/contract 
the bump [17]. Our analysis will emphasize the region of parameter space near where 
bumps are marginally stable to even perturbations. 

2.1. Existence and stability of bumps. We now briefly review existence and 
stability results for stationary bump solutions to the neural field equation (1.1). These 
results are analogous to those presented in [1,35,52]. For transparency, we focus on 
the case of a Heaviside firing rate function (1.3). This allows us to cast bump stability 
in terms of a finite dimensional set of equations, focusing on the evolution of the two 
edge interfaces of the bump [1, 15]. Assuming a stationary solution u{x,t) = U{x), 
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we find (1.1) requires 


U{x)= ( w{x-y)H{U{y)-e)Ay. (2.1) 

Jn 

Given a unimodal bump solution U(x), without loss of generality, we can fix the center 
and peak of the bump to be at the origin a; = 0. In the case of even-symmetric bumps 
U{x) = U{—x) [1], we will have the conditions for the bump half-width a: U{x) > 9 
for X G (—a,a), U{x) < 9 ioi x € fl\[—a,a], and [/(±a) = 9. In this case, (2.1) 
becomes 

/ a px-\-a nx-\-a nx — a 

w{x-y)Ay= / w{y)Ay = / w{y)Ay - / w{y)Ay. 

-a J x—a 0 *^0 

By utilizing the integral function (1.4), we can write the even-symmetric solution 

U{x) = W{x + a)-Wix-a). (2.2) 

To determine the half-width a, we require the threshold conditions U{zta) = 9 oi the 
solution (2.2) to yield 


U(a) = W{2a)= / w{y)dy = 9. 

Jo 

Note that when 9 < Wmax = nia.XxW{x), there will be a stable and unstable bump 
solution to (1.1). When 9 = 9c = Wmax, there is a single marginally stable bump 
solution Uc{x) to (1.1), as illustrated in Fig. 2. IB. Differentiating W{2a) by its ar¬ 
gument yields W'{2ac) = w(2ac) = 0 as an implicit equation for the half-width Uc at 
this criticality. Utilizing the notation of Amari condition (i), we have that Oc = xqI2. 
Note, the relation w(2ac) = 0 is explicitly solvable for Oc for several typical lateral in¬ 
hibitory type weight functions. For instance, in the case of the difference of Gaussians 
w{x) = on x G (— 00 , 00 ) [1], we have Oc = ayJ\n{l/A)/ \2\J— l] 

and 9c = ^ [erf(2ac) — ACTerf(2ac/cr)]. For the “wizard hat” w(x) = (1 — on 

X G (— 00 , 00 ) [12], we have Oc = 1/2 and 9c = e~^. For a cosine weight w(x) = cos(x) 
on the periodic domain x G [—7r,7r] [35], we have Uc = 7r/4 and 9c = 1. 

To characterize the stability of bump solutions to (1.1), we will study the evolution 
of small smooth perturbations eip{x,t) (e <C 1) to stationary bumps U{x) by utilizing 
the Taylor expansion u{x,t) = U{x) + e'ip{x,t) -I- 0(e'^). By plugging this expansion 
into (1.1) and truncating to 0{e), we can derive an equation whose solutions constitute 
the family of eigenfunctions associated with the linearization of (1.1) about the bump 
solution U{x). We begin by truncating (1.1) to 0{e) assuming u is given by the above 
expansion and that the nonlinearity f(u) is given by the Heaviside function (1.3), so 

= -ipix, *) + “ y)H\U{y) - 9)ip{y, t)dy, (2.3) 

and we can differentiate the Heaviside function, in the sense of distributions, by noting 
H{U{x) — 9) = H{x + a) — H{x — a), so 

6{x + a) — 5{x — a) = • ^ - - y- ^ = H'{U{x) — 9)U'{x), 

dx 
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which we can rearrange to find 

H'[U{x) -9)= + (<5(0. + a) + 6{x - a)). (2.4) 

Upon applying the identity (2.4) to (2.3), we have 

^ ^ ^) + 7 [^(a^ + o,)'4’{-<^y t) + w{x - a)iA(a, t)] , (2.5) 

where 7 “^ = |U'(a)| = ■u;(0) — w{2a). One class of solutions, such that ^^(±a,t) = 
■0(410,0) = 0, lies in the essential spectrum of the linear operator that defines (2.5). 
In this case, ip{x,t) = ip(x,0)e~*, so perturbations of this type do not contribute 
to any instabilities of the stationary bump U{x) [24]. Assuming separable solutions 
'ip{x,t) = b{t)ip{x), we can characterize the remaining solutions to (2.5). In this case, 
b'{t) = Xb{t), so b{t) = where A S K, and 

(A + l)'0(a;) = 7 [w{x + a)ilj{—a) + w{x — 0 )^( 0 )]. (2.6) 

Solutions to (2.6) that do not satisfy the condition ip{±a) = 0 can be separated into 
two classes: (i) odd ip{a) = —ipi^—a) and (ii) even ipia) = ’4>{—a). This is due to the 
fact that the equation (2.6) implies the function tp{x) is fully specified by its values 
at a: = ± 0 . Thus, we need only concern ourselves with these two points, yielding the 
two-dimensional linear system 


(A -I- I)0(—a) = 7 [w(O)'0(—a) -I- w(2a)0(a)] 
(A -I- 1)0(0) = 7 [w(2a)0(— o) -I- w(O)0(a)]. 

For odd solutions 0(a) = —0(—o), the eigenvalue 


A 


O 


— 1 -I- 7 [w(0) — w(2a)] = —1 -I- 


u;(0) - w;(2a) _ 
u;(0) - w(2a) ’ 


(2.7a) 

(2.7b) 


reflecting the fact that (1.1) is translationally symmetric, so bumps are marginally 
stable to perturbations that translate their position. Even solutions 0(a) = 0(—a) 
have associated eigenvalue 


Ae 


— 1 -I- 7 [w(0) -I- w(2a)] = —1 -f 


w{0) + w(2a) 
w(0) — w(2a) 


2w(2a) 

w(0) — w(2a)' 


Thus, when 0 < 0^-, the wide bump Og > will be linearly stable to expand¬ 
ing/contracting perturbations since w{2as) < 0 due to Amari’s condition (ii) [1]. The 
narrow bump a„ < Gc is linearly unstable to such perturbations since w{2au) > 0 
due to condition (i). When 0 = 0^, we have w{2ac) = 0 so that Ae = 0 and 
|C/'(±ae)| =u;(0). 

In anticipation of our derivations of amplitude equations, we define the eigenfunc¬ 
tions at the criticality 0 — 0c- Utilizing the fact that |17'(±ac)| = w(0) and the linear 
system (2.7a), we have that the odd eigenfunction at the bifurcation is 


0o(x) 


i(;(0) 


[w(a; - Oc) 


w{x -I- ac)], 


( 2 . 8 ) 
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1 


(2.9) 


and the even eigenfunction is 




;( 0 ) 


){x 


Note, this specihes that V’e(±a) = 
it useful to compute the derivatives 


Oc) + w{x + Oc)] . 

—tjjQ^—a) = 1. Furthermore, we will find 


■^'oix) = [w'{x -a)- w'{x + a)], 

which is even ('0o(—flc) = V'o(®c))) and 

ip'Ax) = [w'{x -a)+ w'{x + a)], 


which is odd {'ip'^{—ac) = —'4’'^{ac)). Lastly, we note that we will utilize the fact 
that, for even symmetric functions, w'(0) = 0, so ip'g{±ac) = '0e(nc) = ~V'e(~nc) = 
w' (2ac)/w{0). 

2.2. Saddle-node bifurcation of bumps. Motivated by the above linear sta¬ 
bility analysis, we now carry out a nonlinear analysis in the vicinity of the saddle-node 
bifurcation from which the stable and unstable branches of stationary bumps emanate. 
Specifically, we will perform a perturbation expansion about the bump solution Uc(x) 
at the critical threshold value 0c- We therefore define 0 — + e 1, so that is 

a bifurcation parameter determining the distance of 9 from the saddle-node bifurca¬ 
tion point. As demonstrated above, the linear stability problem for Uc{x) reveals two 
zero eigenvalues Aq = Ag = 0 associated with the odd and even ipe eigenfunctions 
(2.8) and (2.9), respectively. Our analysis employs the ansatz: 

u{x, t) = Uc{x) + eAe{T)'lpeix) + Ao{t)^o{.x) + £^U 2 {x, t) + 0{e^), (2.10) 


where r = et is a temporal rescaling that reflects the vicinity of the system to a 
saddle-node bifurcation associated with the even expanding/contracting eigenmode 
ijje [51]. Similar expansions have been utilized in the analysis of bifurcations for 
spatial patterns in reaction-diffusion systems [3,49] and neural held models [5,28,53]. 
Upon plugging (2.10) into (1.1) and expanding in orders of e, we hnd that at 0(1), 
we simply have the stationary bump equation (2.1) at 0 = 9c- Proceeding to 0{e), 
we hnd 


0 = Ae(r) 



y)H'{Uciy) - 9c)'ipeiy)dy - tpe{x) > 


so we can use (2.4) to write 


0 = Ac{t) 


1 

w(0) 


(w{x -I- a)'i/>e(—a) -I- w{x 


a)tpe{a)) - i’eix) 


( 2 . 11 ) 


The right hand side of (2.11) vanishes due to the formula for the even (2.9) eigen¬ 
function associated with the stability of the bump Uc(x). At O(e^), we obtain an 
equation for higher order term U 2 '- 


C [Ao^o + -^ 2 ] =A’ci>e + A'c-iIjo + fJ. [ w{x - y)H'{Uc{y) - 0c)dy (2-12) 

Jq 

f w[x-y)H”{Uc{y)-9c)i}e(yfdy, 

^ JQ 
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where £ is the non-self-adjoint linear operator 


Cu{x) = -u{x) -I- / w{x - y)H'{Uciy) - Oc)u{y)dy. (2-13) 

Jo, 

Both ijjo{x) and tpe{x) lie in the nullspace N{C), as demonstrated in the previous 
section by identifying solutions to (2.3). Thus, the tpo terms on the left hand side of 
(2.12) vanish. We can ensure a bounded solution to (2.12) exists by requiring that the 
right hand side be orthogonal to all elements of the nullspace of the adjoint operator 
£*. The adjoint is defined with respect to the inner product 

{Cu,v) = j [£u{x)]v{x)dx = ( u{x)[C*v{x)]dx = {u,C*v). (2.14) 

Jq Jq 

Thus, we find 

C*v{x) = -v{x) + H'{Uc{x) - Oc) [ w{x-y)v{y)dy, (2.15) 

Jn 

defined in the sense of distributions under the inner product given in (2.14). It is 
straightforward to show that (po ■= H'{Uc — dc)fpo and (pf, := H'iUc — 0c)'4’e lia in the 
nullspace of £*. Components of N{C*) are defined by the equation 

v{x) = H'{Uc{x) - Oc) [ w{x - y)viy)dy. (2.16) 

JQ 

To show Po, Pe S we simply plug these formulas into (2.16) to find 


H'{Ucix) - 9c)'ipj{x) = H\Ucix) - 9c) [ w{x-y)H'{Uc{y)-9c)ijjj{y)dy, 

Jq 

for j = o, e, which is true due to the fact that ipo and ipe lie in JV{C). Thus, we will 
impose solvability of (2.12) by taking the inner product of both sides of the equation 
with respect to Po ■= H'{Uc — 0c)V'o and Pc '■= H'{Uc — 9c)'4>e yielding 

0 = (pj,A'oiJc + Ki’o + pw * H'{Uc - 9c) - * [H'\Uc - 9c)^l] ^ , (2.17) 


for j = o, e, where we have defined the convolution w*F = w(x — y)F(y)dy. Due to 
odd-symmetry, terms of the form {FI'[Uc — 9c)’4>j, ’ipk), j 7^ k, vanish. In a similar way, 
the term {H'{Uc — 9c)'4>o,w* H'{Uc — 9c)) vanishes due to odd-symmetry. Isolating the 
temporal derivatives Aj in (2.17), we find that the amplitudes Aj (j = o, e) satisfy 
the following fast-slow system of nonlinear differential equations 


dAo {po, w * [H"{Uc - 9c)^l] ) 2 

dt 2{po,^o) ’ 

dA _ _ {pc,w*H'{Uc-9c)) {pe,w-^[H"{Uc-9c)i)l]) 2 

dr ^ {Pe,1pe) 2{pc,i!e) 


(2.18a) 

(2.18b) 


With the system (2.18) in hand, we can determine the long term dynamics of 
the amplitudes as the bifurcation parameter p is varied. We begin by computing the 
constituent components of the right hand sides, using properties of the eigenfunctions 
ipo and ipc- To start, we will compute the second derivative F["{Uc—9c), which appears 
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in the coefficient of the quadratic term A\. Differentiating the function H{Uc{x) — 6c) 
twice with respect to a;, using the chain and product rule, we find the following formula 

d^HiUcix) - 9c) 


dx"^ 


= {U',{x))^H"{Uc{x) - dc) + U'Jix)H'{Ucix) - 6c) 


= {U'{x))^H"iUcix) - 6c) + 


C/"(x) dH{Uc{x) - dc) 


U'cix) 


dx 


where we have applied the identity (2.4) for the first derivative H'{U—6). Rearranging 
terms, we find that 

1 d^H{Uc{x)-6c) U':{x) , , ,, 

H {Uc ^^2 |C/'(a)|3 + 

We can further specify the formula (2.19) by differentiating = ^{x + Uc) — 

6{x — Uc) with respect to x to yield 

d^H(Uc - 6c) 


dx^ 


= S'(x + Uc) — s' (x — Uc), 


where S'(x — xo) is defined, in the sense of distributions, for any smooth function F(x) 
by using integration-by-parts [32] : 

I S'{x — xo)F{x)dx = — j S{x — Xo)F'{x)dx =—F'{xo). 

Ja JQ 

Furthermore, we note that the spatial derivatives |17'(±ac)| = rc(0) and U”{x) = 
w'{x + Uc) — w'{x — ttc). Even symmetry of w(x) mandates that w'(x) = —w'(—x) 
and w'(0) = 0, so 17"(±ac) = w'{2ac). Thus, we can at last write 

cjff/rr ^ ^ _ ^'(a; + ac)-^'(a;-ac) w'{2ac)[S{x + Gc) + S{x - Uc)] 

\^c ^c) — (c\\2 ^r\\3 ' 

w[0)^ w{\jy 

Computing the inner products in (2.18) then simply amounts to evaluating the inte¬ 
grals in the sense of distributions. First, we use (2.4) to note 

= I ^j{x)^H'{Ucix) - 6c)dx = 7 [V’i(ac)^ -b V't(-ac)^] = 

for j = o, e, since -ipeidzac) = ipoiac) = —V'o(—Oc) = 1- Furthermore, 

{(Pc,w * H'{Uc-dc)) = [ [ w{x - y)ipcix)H'{Ucix) - 6c)H'{Uciy) - 6c)dydx 
Jfl JQ 


= 7 


w{x + a) 


.a=iac 


Ipeix) 


S(x + a) 
2 


a=±a, 


da; 


= 7^ [ipeiac) + ipei-ac)] ■ [rc(0) -b w(2ac)] = 


.( 0 )’ 


( 2 . 21 ) 


where we have utilized tpe{d^O‘c) = 1 and w(2ac) = 0. Finally, we compute the 
quadratic terms using the identity (2.20), starting with 

{ipo,w* [H"{Uc - 6c)'ipe]) = [ f “ y)'fo{x)H"{Uc{y) - 6c)ipe{y)'^dydx 

JQ Jn 

= 7 X! / ''xia-y)H"{Uciy)-6c)tpe{y)'^dy, 

a—ihac ^ 


( 2 . 22 ) 
















and we note that individual terms under the integral from the sum dehning (2.20) are 

/ w{-ac - y)S'{y + ac)ipe(,y)'^‘iy = ■u;'(0)'(/'e(-ac)^ - 2ui(0)'!/'e(-ac)V'e{-ac) 

Jn 

= 2w'(2ac), 

/ w{ac - y)S'{y + ac)ipe{y fdy = w'(2ac)V'e(-ac)^ - 2i(;(2ac)V'U-ac)V'e(-ac) 

Jn 

= w'{2ac), 

/ w{-ac - y)6'{y - ac)'ipe{y)'^dy = w'(-2ac)'i/'e(ac)^ - 2w{2ac)tp'e{ac)tpe{ac) 

Jn 

= -w'(2ac), 

/ w(ac - y)^'(y - ac)'tpe(yfdy = w'(0)'ipe(acf - 2w(0)V'e(ac)V'e(ac) 

Jn 

= —2w'{2ac), 

for the terms involving the distributional derivative 6'(x — Xq), whereas the terms 
involving S(x — xq) are 

/ w(-ac - y)5{y + ac)'il^e{yfdy = 'u;(0)V'e(-ac)^ = w(0), 

Jn 

/ w{ac - y)6{y + ac)tj}e{y)'^dy = w{2ac)'ipe{-ac)'^ = 0, 

Jn 

/ w{-ac - y)6{y - ac)'(/'e(2/)^d?/ = w(2ac)'ipe{acf = 0, 

Jn 

/ wiac - y)S{y - ac)'il}eiyfdy = ui(0)V’e(ac)^ = w(0). 

Jn 

Thus, each integral term 

w{-ac - y)H"{Uc{y) - Oc)tpe{yfdy 

f w{ac - y)H"{Uc{y) - 0^)tpe{yfdy = ■ 

Jn w{oy 

Finally, using the fact that il’oio.) = ~4’o{~o) = 1, we find that the two terms in the 
sum of (2.22) cancel and the integral vanishes. Thus, {(fio,w * [H''iUc — 0c)V’e]) = 0; 
so Ao{t) = Ao is constant. On the other hand, computing the quadratic coefficient in 
the equation for A^, we have 

{(fe.w* [H"{Uc - 0a)'ipl]) = [ [ w{x - y)ipeix)H"{Uciy) - 0c)i^e{y)‘^dydx 

J Q, J n 

= 7 X! f w{a-y)H"{Uc{y) - 0c)ipe{y)'^dy. 

a—dzac ^ 

(2.25) 


(2.23) 

(2.24) 


The integrals in (2.25) are identical to those in (2.22), so it is straightforward to 
compute, using (2.23) and (2.24) that 


{(Pe,W* [H"{Uc - 6»c)'!/'e] ) = 7 


2w'{2ac) ^ 2w'{2ac) 


,{Qf 


,{Qf 


Aw'{2ac) 

w{QY 
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(2.26a) 


Thus, we can at last compute all the terms in (2.18), specifying that 


dt 

dyl, 

dr 


= 0 , 


= -M- 


|w'(2ae)| 

u;(0)^ 


Aeirf 


(2.26b) 


where we have noted the fact that w'{2ac) < 0 due to Amari’s conditions (iii) and 
(iv) on the weight function w(x) [1]. 

Equation (2.26a) reflects the translational symmetry of the original neural field 
equation (1.1), so bumps are neutrally stable to translating perturbations ipo regard¬ 
less of the bifurcation parameter /r. On the other hand, as the bifurcation parameter 
/r is changed, the dynamics of the even eigenmode ipe reflect the relative distance to 
the saddle-node bifurcation at which point bumps are marginally stable to expand¬ 
ing/contracting perturbations. When /r < 0, there are two fixed points of equation 
(2.26b) at Ae = ±w{0)^/\fi/w'{2ac)\, corresponding to the pair of emerging station¬ 
ary bump solutions which are wider (-I-) and narrower (—) than the critical bump 
Uc- As expected, the wide bump is linearly stable since a linearization of (2.26b) 
yields A+ = ■ u>'(2ac)|/w(0) < 0, and the narrow bump is linearly unstable 

since A_ = +V\i^ ■ w'{2ac)\lw(0) > 0. Crossing through the subcritical saddle-node 
bifurcation, we find that for /i = 0, there is a single fixed point Ag = 0, which is 
marginally stable, since Aq = 0. 

Lastly, note when fj, > 0, there are no fixed points of the differential equa¬ 
tion (2.26b). However, starting at the initial condition Ae(0) = 0 (correspondingly 
m(x, 0) = Uc{x)), we find that the dynamics of the amplitude Ae(T) are strongly de¬ 
termined by the ghost of the fixed point at Ag = 0 [51]. Note in Fig. 2.2A that the 
transient bump retains a shape much like that of the critical bump for an appreciable 
period of time before extinguishing. Trajectories of the full system (1.1) evolve more 
slowly when the distance to the bifurcation \6 — 0g| = |/r|£^ is smaller. Solving for 
Ag(T) in this specific case and reverting the the original time coordinate t = r/e, we 
hnd 

= ~ (eV/r- |w'(2ag)|V'f^(0)) . (2.27) 

^y\w'(2ac)\ ^ ^ 


Thus, the residence time tb in the bottleneck, or neighborhood of the ghost of the 
fixed point Ag = 0, is given by the amount of time it takes for Ag(t) to traverse 
to some set value. Of course, this is dependent on the bifurcation parameter g,. For 
illustration, we examine how long it takes until Af.(tb) = — 1. By explicitly focusing on 
the region where |Ae(tb)| < 1, we are roughly restricting to the time interval during 
which \u{x,t) — Uc{x)\ = 0{e), where we would expect the expansion (2.10) to be 
valid. Using the formula (2.27), it is straightforward to find that 


e-y/pA]^'7(2^ \ j 


(2.28) 


We compare this formula to the results of numerical simulations in Fig. 2.2B, utilizing 
the difference of Gaussians weight function w{x) = e““ — Ae“^ on a: € (—oo, oo). 
Comparisons are made by noting that when Ae(tb) = —1, then u{x,t) « Uc{x) — 
eipeix), so that the peak of the activity profile will be 

m(0, tb) « C/g(0) - eV'g(O) = W(ag) - W(-ag) - 
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Fig. 2.2. Slow passage of bumps on x G (— 00 , 00 ) when w(x) = e~^ — Ae~^ . (K) Slow 
passage of a transient bump by the ghost of the critical solution Uc{x) when 9 = 6c for e = 0.1 
(fj, = 1), A = 0.4, and a = 2. (BJ The peak of the bump u{0,t) slowly decreases in amplitude until 
breaking down quickly in the vicinity of Ae{t) = — 1. Note the theoretical formula for the amplitude 
(solid line) given by (2.27) matches the numerical simulation (dashed line) in the slow passage 
region. (Q) Amplitude of the even mode Ae(t) slowly decreases with time. The duration of the 
bottleneck increases as the distance to the bifurcation is decreased by reducing s. (Q) Comparison 
of the theory (solid) given by (2.28) to the numerically computed (dots) duration in the bottleneck 
(the crossing Ae(ti,) = —Ij. 




Fig. 2.3. Slow passage of a bump on x £ [— 7r,7r] for a cosine weight function w(x) = cos{x). 

(A) Amplitude of the even mode Ae{t) slowly decreases with time. Numerical simulations (dashed 
lines) of (1.1) are compared to the trajectory \/2{l — £:tan(£t)) determined by theory (solid lines). 

(B ) The duration of bottleneck increases as the distance to the bifurcation is decreased. Simulations 
(dots) are well fit by the theory = '7r/[4£:] (solid line). 
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Notice in Fig. 2.2C,D that, as predicted, the time spent in the bottleneck increases 
as the amplitude of the small parameter e is decreased. The attracting impact of 
the ghost is stronger when the parameters of the system lie closer to the bottleneck. 
For further comparison, we consider the case w(x) = cos(x) in Fig. 2.3. In this case 
the constituent functions Oc = 7r/4, i(;(0) = 1, and w{2ac) = —1. Furthermore, by 
setting ^ = 1 the formulas for the amplitude (2.27) and residence time (2.28) simplify 
considerably to Ae{t) = — tan(£t) and 4 = 7r/[4£]. 


2.3. Amplitude equations for smooth nonlinearities. Our nonlinear anal¬ 
ysis in the case of Heaviside nonlinearities f{u) = H{u — 9) made extensive use of the 
specific form of the distributional derivatives. Inner products with these functions 
lead to dynamical equations focused on a finite number of discrete points in space, 
rather than over the spatial continuum x S H. Here, we show it is straightforward 
to extend this analysis to the case of arbitrary smooth nonlinearities f(u). There 
are several detailed analyses of stationary bumps in neural field with smooth firing 
rate, showing a similar bifurcation structure to that presented in Fig. 2.1: a stable 
and an unstable branch of bump solutions annihilate in a saddle-node bifurcation as 
the threshold of the firing rate function is increased. We refrain from such a detailed 
analysis here and refer the reader to these works [13,18,35,37,41,52]. Again, defining 
9 = 9c + /is^, e 1, so /r determines the distance of 9 from the bifurcation and on 
which side of 9c it lies. Following our previous analysis, we utilize the ansatz (2.10) 
and rescale time t = ex. In this case, ipoix) and ipeix) will still be odd and even 
eigenmodes associated with the linear stability of stationary bump solutions to (1.1). 
At the criticality 9 = 9c, their associated eigenvalues will be Ao = Ae = 0, as in the 
case of Heaviside firing rates [52]. Expanding (1.1) in orders of e using the ansatz 
(2.10) yields a similar amplitude equation to (2.18) at 0{e^). Again, we apply solv¬ 
ability conditions to the equation for U 2 - After canceling odd terms and isolating the 
derivatives A', we find the amplitudes Aj satisfy the system: 


dAo {(Po,w* [f'{Uc)'tp'^]) 2 

dt 2{^c,i’o) ’ 

dAe ^ _ {ipc,W* f{Uc)) {ifc,W*[f''{Uc)i^l]) 2 

dr ^ {^e,'Pe) 2{(pc,^e) 


(2.29a) 

(2.29b) 


We can derive the coefficients in the system (2.29) by computing the inner prod¬ 
ucts therein. To do so, we must choose a specific nonlinearity, such as the sigmoid 
(1.2), and a weight kernel. For illustration, we consider the cosine kernel w(x) = cos(x) 
on the ring x S H = [—tt, tt] with periodic boundaries. As shown in previous stud¬ 
ies, the bump solution Uc{x) = Ac cos x while the eigenmodes ipo{x) = sin(x) and 
tpeix) = cos(x) [26,35,52]. Since £tpj = 0 for j = o,e, this means 

/ TT pTT 

cos{x - y)f {Ac cos{y))sm{y)dy = sinX / sin^(y) f {Ac cos y)dy, 

-TT J —TT 


where we have used cos(x — y) = cos x cos y + sin x sin y, and 

/ TT pTT 

cos{x - y)f {Ac cos{y)) cos{y)dy = cos X / cos^(2/)/'(Ac C0S2/)d2/, 

-TT J —TT 


so that we can write 


/ TT pTT 

sin^ {y)f'{Ac cos y)dy = 1, / cos^{y)f{AcCOsy)dy = 1. 

-TT J —TT 


(2.30) 
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The identities (2.30) allow us to compute 

{tpo, ipo) = / f'{Ac cos{y)) sin(?/)^dj/ = 1, 

J —TT 

and 

(V5e,V'e)= / f'{AcSui{y))cos{yf‘dy = l. 

J —TT 

Furthermore, 

{ipo,w*[f'{Uc)Tpl]) = f f"{Uciy))'ipe{yf f cos{x - y)f'{Ac cos{y))sm{y)dxdy 

J —TT J —TT 

= / /"(17c(2/))cos(?/)^sin(j/)dy = 0, (2.31) 


where the last equality holds due to the integrand being odd. Thus, the equation 
(2.29a) reduces to A'^{t) = 0, so Ao{t) = Ag. Now, we can calculate the coefficients 
of the Ac amplitude equation. First by utilizing the fact that w{x — y)(pe{y)d-y = 
ipeix), we can compute 

{(Pe,w * f{Uc)) = f f {Accos{x))cos{x)dx = {ipeA)- {2.32) 

J —TT 


Lastly, we can simplify the integrals in the quadratic term by again making use of the 
identity w{x - y)ipe{y)dy = 'tpe{x), so 


{ificW* [f'{Uc)tpl]) = j f'{Ac cos{x))cos^{x)dx = {f'{Uc),tpl). 
J —TT 


(2.33) 


so we can simplify (2.29b) to 

d A 1 

^ = -y{y,cA) + -{r{Uc),^Pl)Ac{T)\ (2.34) 

3. Stochastic neural fields near the saddle-node. We now study the impact 
of stochastic forcing near the saddle-node bifurcation of bumps. Our analysis utilizes 
the spatially extended Langevin equation with additive noise (1.5). Guided by our 
analysis of the deterministic system (1.1), we will utilize an expansion in the small 
parameter £, which determines the distance of the system from the saddle-node. To 
formally derive stochastic amplitude equations, we must specify the scaling of the 
noise amplitude e as it relates to the small parameter e, as this will determine the 
level of the perturbation hierarchy wherein the noise term dW will appear. We opt for 
the scaling e = as this introduces a nontrivial interaction between the nonlinear 
amplitude equation for Ac and the noise. 

It is important to note that our derivations are only carried up to O(e^) in the 
hierarchy of the regular perturbation expansion in e. Were we to continue this expan¬ 
sion further, we would likely find that the e = noise term does indeed shift the 
location of the bifurcation at higher order as in [2,30]. Thus, as the amplitude of noise 
is increased, the validity of the expansion we derive here will begin to break down, 
since the terms beyond 0{e^) will have a more substantial effect on the dynamics. 
Hence, the results we derive in this section are valid for small noise levels only. An 
understanding of the effects of larger noise terms, employing scalings e = with 
p < 5/2, warrants further study which is beyond the scope of our current work. 
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3.1. Stochastic amplitude equation for bumps. Motivated by our quanti¬ 
tative analysis in the noise-free case, we rescale time in the stochastic term of (1.5) 
using T = et, so 


du{x, t) 


—u(x, t) 



y)f{u{y,t))Ay 


dt + e^dW{x, r), 


(3.1) 


where dW{x,T) := s/edW{x, t) is a rescaled version of the Wiener process dW 
that is independent of £ [22]. We then apply the ansatz (2.10) once again and take 
Heaviside firing rate functions (1.3), thus finding (2.11) at 0{e). The 0(e) equation 
is satisfied due to the fact that V'e € A/’(£), where £ is the linear operator given by 
(2.13). Finally, proceeding to 0{e^), we find 


£ [Ao^po + U 2 ] dt =dAeipe + dAoipo + M / w{x - y)H'[Uc{y) - 0c)dydt (3.2) 

Jo, 

- 4^ / w(x-y)H”{Uc{y)-6c)ipe{y)'^dydt + dW. 

As before, the ^|Jo terms on the left vanish since Cipo = 0, and we ensure a bounded 
solution to (3.2) exists by requiring the inhomogeneous part is orthogonal to ipo, ‘fe G 
A/'(£*), where £* is the adjoint linear operator given by (2.15). Taking inner products 
yields 

0 = {(pj,dAe{T)p)e{x) + dAo(t)'0o(a;) + yw * H'{Uc - 0c)dt (3.3) 

* [H"{U, - 0e)V^e] di + dw^ , 


for j = o, e. Isolating temporal derivatives, we find the amplitudes Ao{t) and Aeij) 
obey the following pair of nonlinear stochastic differential equations 




2(^0; 


(^cV'o) 


dAe(r) = - /i 




{‘Pe,1pe) 


2(^65 V^e) 


Aeir)^ 


(3.4a) 

(3.4b) 


_ (v?e,dW) 

(£e,V'e) ’ 

Utilizing the formulas for H'{Uc — dc) (2.4) and H"{Uc — Oc) (2.19) we derived in the 
previous section, we can simplify the expressions in (3.4). Additionally, we make use 
of the fact that 


dyVo(T) 


{iPo,dW) 

(£o,V'o) 


d>Ve(T) 


(tpe,dlU) 

(£e,V’e) 


'ipo(-ac)dW(-ttc, t) + '0o(ac)dlU(oc, t) 


2 L 


dlU(-flc, r) - dW (oc, r) 


V'e(-ac)dlU(-ac,r) +'0e(ac)dlU(ac,r) 
dlU (oc, r) -I- dW (-Oc, r) 
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Fig. 3.1. Noise-induced extinction of bumps in the stochastic neural field (1.5) on x £ [—7r,7r] 
for a cosine weight w(x) = cos(a;). (A.) A single realization of the equation (1.5) with the initial 
condition w(a:,0) = Uc{x) = y/2cos{x) leads to a stochastically wandering bump that eventually 
crosses a separatrix at t ^ 70, leading to extinction. The noise-free system possesses a stable 
bump solution since fi = —0.2 < 0; e = 0.4. fB) The large deviation can easily be detected by 
tracking max^ii(a:;, t), which departs the bottleneck of the noise-free system, whose lower bound lies 
at max^ [Uc{x) — sTpeiO)] = \/2(l — e). 


utilizing the fact that (dlU(a:, r)dlU(j/, r')) = C{x — y)S{T — T')dTdT', it is straightfor¬ 
ward to compute the variances (>Vo('r)^) = DoT = ((7(0)—C'(2ac))'r/2 and (>Ve('r)^) = 
DeT = ((7(0) + C{2ac))T/2. Clearly, for spatially flat correlation functions C{x) = (7, 
noise will have no impact on the odd amplitude Ao{t) since Do = 0. Thus, (3.4) 
becomes 


dy4o(t) = ^/edWoit), (3.5a) 

dAe{r) = -fidr - + dVVe(T), (3.5b) 

w[0)^ 

where we have converted the noise term in (3.5a) back to the original time coordinate: 
dWo(t) = d'Wo{£t)/y/e [22]. Note that in equation (3.5a), we essentially recover the 
diffusion approximation of the translating mode of the bump {Ao{t)'^) = sDot, which is 
analyzed in [35]. Equation (3.5b) is a stochastic amplitude equation, so that the noise 
term dWe is projected onto the direction of the neutrally stable even perturbation 'ipo- 

3.2. Metastability and bump extinction. To analyze the one-dimensional 
nonlinear SDE (3.5b), we further rescale the equation by setting A := 

d^(T) = — [m-I-A( t)^] dt-I-dyV(r), (3.6) 

where m := Thus, the effective diffusion coefficient of the rescaled noise 

term is (VV('r)^) = Dt = w'{2acY{C{Q) + (7(2ac))T/ [2w(0)^]. Note the rescaled 
equation (3.6) has an effective potential [42,51]: 

V{A) = ^+mA, (3.7) 

the derivative V'{A) of which yields the deterministic part of the right hand side. As 
the bifurcation parameter m is varied, the potential exhibits a minimum (at A = .Jm) 
and a maximum (at A = —yAri) when m < 0, a saddle point (at A = 0) when m = 0, 
and no extrema for m > 0 (Fig. 3.2A). For all parameter values m, the state of the 
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Fig. 3.2. (A) Potential function (3.7) associated with the stochastic amplitude equation (3.6) 
has zero (m > 0); one (m = 0); or two (m < 0^ extrema - associated with equilibria of A = —m — A^. 
When m < 0, crossing the saddle point requires stochastic forcing. ) Mean time tf, until bump 
extinction is approximated by a mean first passage time problem of the stochastic amplitude equation 
(3.6). Numerical simulations (circles) of the full system (1.5) are well approximated by this theory 
(line) given by (3.12) for e = 0.6. 


stochastic system (3.6) will eventually escape to the limit A —> —oo as t — )■ oo. Such 
trajectories were observed in the noise-free system in the case m > 0, as demonstrated 
in Fig. 2.2 of the previous section. However, we show here that noise qualitatively 
alters the dynamics of the system, so its state will not remain in the vicinity of the 
stable attractor (at A = y/m) when m < 0. 

As before, we study the problem of bump extinction using the stochastic am¬ 
plitude equation (3.6) in the case m > 0. We show that the noise decreases the 
average amount of time until an extinction event will occur. For clarity, we as¬ 
sume the initial condition A(0) = 0 (correspondingly u{x,0) = Uc{x)). We take 
the bottleneck to be the region Ae G [—1,1], which in the rescaled variable is A S 
[—|■u;'(2ac)|/u^(0)^, |ic'(2ac)|/rc(0)^]. The residence time u in the bottleneck is given 
by the amount of time it takes for A to escape this region. We can determine the 
statistics of Tf, by considering it as a first passage time problem. 

Let p{A, t) be the probability density for the stochastic process A{t) given the 
initial condition A(0) = Aq. Then the corresponding Fokker-Planck equation is given 

dp ^d[{m + A^)p{A,T)] Dd^p{A,T) _ dJ{A,T) 

dr dA 2 dA^ ~ dA ’ ^ 


where 


J(A,r) =- (m-h A2 )p(A,t), (3.9) 

and p{A, 0) = 5{A — Aq). We focus on the three different scenarios discussed above. 
First, if m < 0, there there is a single stable fixed point of the deterministic equation 
A = —m—A^ at A = .Jm and a single unstable fixed point at A = —y/m. The basin of 
attraction of A = y/m is given by the interval {—y/m,oo). When D > 0, fluctuations 
can induce rare transitions on exponentially long timescales whereby A(t) crosses the 
point A = —y/m, leaving the basin of attraction. For the non-generic case m = 0, the 
timescale of departure scales algebraically [50]. When m > 0, noise simply modulates 
the flows of the deterministic equation A = —m — A^, leading to an average speed-up 
in the departure from the bottleneck. In general, we consider solving the first passage 
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time problem as an escape from the domain {—a, oo) where a := (equivalently 

where Ag = —1) [22]. To do so, we impose an absorbing boundary condition at —a: 
p{—a,T) = 0. Now let T{A) denote the stochastic first passage time for which (3.6) 
first reaches the point —a, given it started at A € (—a, oo). The first passage time 
distribution is related to the survival probability that the system has not yet reached 
—a: 

pOO 

S{t)= / p{A,T)dA, 


which is S{t) := Pr(t > T{A)), so the first passage time density is [22] 


dS 


dp 


F{t) = -— = - —{A,T)dA. 


dr 


dr' 


Substituting for the expression for dp/dr using the Fokker-Planck equation (3.8) and 
the formula for the flux (3.9) shows 


F(r) 



dJ{A, t) 
dA 


d^ = 


-J{-a,T), 


where we have utilized the fact that lim^i-^oo J^A, r) = 0. Thus, the first passage time 
density F{t) can be interpreted as the total probability flux through the absorbing 
boundary at A = —a. To calculate the mean first passage time T{A) := {T{A)), we 
use standard analysis to associate T{A) with the solution of the backward equation 
[ 22 ]: 


-{m + A )— + 


dr D d^T 


dA 2 dA^ 


= -l, 


(3.10) 


with the boundary conditions T{—a) = 0 and T'(oo) = 0. Solving (3.10) yields the 
closed form solution 


r{A) = 



Hz) 

Hy) 


dzdy, 


(3.11) 


where 


(j){A) = exp 


2 [v(-a) - v(A)] 
D 


and V{x) is the potential function (3.7). Explicit expressions for the integral (3.11) 
can be found in some special cases [42,50]. For our purposes, we simply integrate 
(3.11) numerically to generate theoretical relationships between the mean first passage 
time and model parameters. For comparison, we focus on the case the weight function 
w{x) = cos(x) and the correlations C{x) = cos(x), so that Uc{x) = •\/2cos(a::), Oc = f, 
w{0) = 1, w'{2ac) = —1, C'(O) = 1, and C{2ac) = 0. Therefore, a = I, m = p, 
D = 1/2. This allows us to write the formula (3.11) at A = 0 as 


/ O poo 

/ exp 
-1 Jy 




+ piz-y) 


dzdy. 


(3.12) 


Lastly, note that by rescaling time t = er, we have that the mean first passage time 
in units of t will be ib = T(0)/e. We compare our theory (3.12) with the results of 
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numerical simulations of the full stochastic neural field (1.5) in Fig. 3.2B. Note there 
is some discrepancy between our numerical simulations and theory as m is decreased. 
One of the primary reasons for this deviation is likely because of the moderate level 
of noise (e = 0.6) used in comparison to the small parameter assumption (e ^ 1) 
using in the theory we have developed. Any minor mismatch will be exacerbated by 
the fact that mean first passage times for escape problems depend exponentially on 
parameters like noise amplitude and well depth, as in (3.12). Nonetheless, the theory 
does provide a rough estimate of the mean first passage times for smaller values of 
the parameter m. 

4. Discussion. We have developed a weakly nonlinear analysis for saddle-node 
bifurcations of bumps in deterministic and stochastic neural field equations. While 
most of our analysis has focused upon Heaviside firing rate functions, we have also 
demonstrated the techniques can easily be extended to arbitrary smooth nonlineari¬ 
ties. In the vicinity of the saddle-node, the dynamics of bump expansion/contraction 
can be described by a quadratic amplitude equation. For deterministic neural fields, 
this low dimensional approximation can be used to approximate the trajectory and 
lifetime of bumps as they slowly extinguish. To do so, we focused on the initial time 
epoch in the bottleneck surrounding the ghost of the critical bump Uc(x). In stochas¬ 
tic neural fields with appropriate noise scaling, a stochastic amplitude equation for 
the even mode of the bump can be derived. Importantly, we must choose the noise 
amplitude to scale as e = in order for the noise term to appear in the stochastic 
version of the quadratic amplitude equation. We then cast the lifetime of the bump 
in terms of a mean first passage time problem of the reduced system, which is valid 
for the noise scaling we have chosen. 

Our work extends a variety of recent studies that have derived low-dimensional 
nonlinear approximations of neural field pattern dynamics in the vicinity of bifur¬ 
cations [5,7,20,30,35,36]. As in our work, most of these previous studies derived 
approximations where the location of the bifurcation was unaffected by noise terms. 
On the other hand, Hutt et al. showed that noise can in fact shift the position of Tur¬ 
ing bifurcations in neural fields, and the amplitude of the bifurcation threshold shift 
was proportional to the noise variance [30]. Were we to have carried the hierarchy 
out to higher order, we would have found such a shift in the case we studied. Note, it 
was necessary in our work to apply a specific noise scaling (£®/^), as compared to the 
distance from criticality (£^), in order for the noise to simply appear as a modification 
of the even mode amplitude equation. Were we to have selected noise of larger am¬ 
plitude, this could have induced bifurcation shifts at lower order, analogous to that 
found in [30]. Another potential future direction would be to consider the impact of 
axonal propagation delays [29] on the dynamics close to the saddle-node. As demon¬ 
strated in this work, the neural field (l.I) is quite sensitive to small perturbations 
near criticality, so delays may alter the duration of the bottleneck or even shift the 
saddle-node bifurcation point. In previous work [36], we derived amplitude equations 
describing propagation-generating drift bifurcations that arise when linear adaptation 
is incorporated into (I.l). We anticipate that similar analyses might be performed on 
networks with synaptic depression [33], although piecewise smooth methods may be 
necessary in the case of Heaviside firing rates [34]. Lastly, we note there have been 
recent efforts to systematically derive macroscopic descriptions of neural activity from 
recurrently coupled spiking networks [9,39,43]. Such alternative descriptions can also 
capture the form of steady state solutions like stationary bumps [39,46]. Bumps 
in these models do exhibit saddle-node bifurcations similar to those observed in the 
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Amari model [1]. Thus, extending our methods to such models would merely require 
knowledge of the stationary equations defining bump solutions. Perturbation analysis 
along with solvability conditions could then yield the coefficients of the amplitude 
equations near the saddle-node. 
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